arXiv: 1502.01809v3 [hep-lat] 14 Apr 2015 


Preprint number: KYUSHU-HET-149 


Universal formula for the flavor non-singlet 
axial-vector current from the gradient flow 

Tasuku Endo^, Kenji Hieda^, Daiki Miura^, and Hiroshi Suzuki'^’* 

Department of Physics, Kyushu University, 6-10-1 Hakozaki, Higashi-ku, Fukuoka, 812-8581, 
Japan 

* E-mail: hsuzuki@phys.kyushu-u.ac.jp 


By employing the gradient/Wilson flow, we derive a universal formula that expresses 
a correctly normalized flavor non-singlet axial-vector current of quarks. The formula is 
universal in the sense that it holds independently of regularization and especially holds 
with lattice regularization. It is also confirmed that, in the lowest non-trivial order of 
perturbation theory, the triangle diagram containing the formula and two flavor non¬ 
singlet vector currents possesses non-local structure that is compatible with the triangle 
anomaly. 
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1. Introduction 


In this paper, we consider a well-studied problem, how to construct a correctly (or canoni¬ 
cally) normalized flavor non-singlet axial-vector current of quarks 0 in a new light using the 
gradient/Wilson flow MS An axial-vector current j^^{x) is said to be correctly normal¬ 
ized, if it fulhlls the Ward-Takahashi relation associated with the flavor chiral symmetry. 
That partially conserved axial current (PCAC) relation is 


= - x)'y5t^ {'il:{y)'4}{z)) - 6^{z - x) (V’(y)V^(2)) (1.1) 


where V' and ^ are quark and anti-quark helds and (A = 1, ..., Nj — 1) denotes the gen¬ 
erator of the flavor SU{Nf) group. Mq is the (flavor-diagonal) bare quark mass matrix. This 
relation says that the axial-vector current j^^{x) generates the axial part of the flavor sym¬ 
metry SU{N f)L X SU{Nf)ii in the correct magnitude. The construction of such a composite 
operator is, however, not straightforward because a regulator that manifestly preserves the 
chiral symmetry does not easily come to hand. The only known explicit exaimles of such 
chiral-symmetry-preservingregularization are the domain-wall lattice fermiori mMg] and the 
overlap lattice fermion [^,10], both satisfy the Ginsparg-Wilson relation 11- 13|. 

For example, if one uses dimensional regularization with complex dimension L) = 4 — 2e 
and 


75 = 70717273, 


( 1 . 2 ) 


a simple expression for the axial-vector current 




(1.3) 


does not fulfill Eq. Instead one hnds (see Ref. [l^ . for examp lej^ 


{df, [i/(x)7/,75t^i/(x)] ijj{y)ijj{z)) - {'ijj{x)'y5{t^, Mo}'ijj{x)' 4 >{y)'ijj{z)) 

= - x)'y5t^ (i/(y)i/(2;)) - S^(z - x) (i/(y)i/(2;)) 75^"^ 

+ ^£^<^ 2(^)4 (<9;, [Vi(x)7/,75t^'0(a:)] ^p{y)^p{z)) 

- ■(£^^2(R)4 (?/(x)75{t^, Mo}i/(x)i/(y)i/(z)) -h 0{g^) (1.4) 


^ See Ref. [J and references cited therein for various methods. 

2 A strategy to determine a non-perturbative renormalization constant of the axial-vector currmt 
on the basis oQhe axial Ward-Takahashi relation in the flowed system has been developed in Ref. [J] . 
See also Ref. [5| for a detailed study of the axial Ward-Takahashi relation in the flowed system. In 
these papers, the flowed fields are employed as a “probe” rather than to construct the axial-vector 
current itself. In the present paper, we instead construct the axial-vector current from the flowed 
fields through the small flow-time expansion. 

^ Here, go is the bare gauge coupling. This expression is for a fermion in a generic gauge rep¬ 
resentation i? of a gauge group G. The quadratic Casimir C 2 {R) is defined by = —C' 2 (i?)l 

from anti-Hermitian generators T“ of G; we normalize generators as trfl(T“T^) = —T{R)6°‘^. For 
the fundamental representation A of G = SU{N), the conventional normalization is T{N) = 1/2 
and G 2 (A) = (A2 - 1)/(2A); C 2 {R) = 4/3 for quarks. 


2 









for e —)■ 0. This relation shows that under dimensional regularization, the correctly normal¬ 
ized axial-vector current is not Eq. (11.31) but 


= 


l-^C2{R)^ + 0{gt) 




in conjunction with a redehnition of the pseudo-scalar density, 

-,2 


l-j^C2{R)^ + 0{gt) 




(1.5) 

( 1 . 6 ) 


The relation between the correctly normalized axial-vector current and a bare axial-vector 
current is regularization-dependent and generally receives radiative corrections in all orders 
of perturbation theory. If one changes regularization (to lattice regularization with a particu¬ 
lar discretization of the Dirac operator, for example), one has to compute the corresponding 
relation anew. 

In the present paper, instead, we derive a single “universal formula” that is supposed to 
hold for any regularization. Our result is 



1 + 


(47r)2 


C2{R) 


-i+ln(432) 


I x{t,x)jf,j5t^x{t,x). 


(1.7) 


Here, g{q) is the runniM gauge coupling in the minimal subtraction (MS) scheme at the 
renormalization scale gjj x(T x) and x(t, x) are quark and anti-quark fields evolved W flow 
equations which will be elucidated below. Because of theorems proven in Refs. [i, 3] , the 
composite operator in the right-hand side of Eq. (ini) is a renormalized quantity although 
it is constructed from bare quark fields in a well-defined manner. As far as one carries out 
the parameter renormalization properly, such a renormalized quantity must be independent 
of the chosen regularization. In this sense, the formula is “universal.” The formula is also 
usable with lattice regularization and could be applied for the computation of, say, the pion 
decay constant. 

Since g{l/V^) —)• 0 for t —)• 0 because of the asymptotic freedom, the formula (11.71) can 
further be simplified to j^^{x) = limt_>.o x(^5 x). Although this is mathematically 

correct, practically one cannot simply take the t —)• 0 limit in lattice Monte Carlo simulations 
for example. In Eq. (ILZD, it is supposed that the regulator is removed (after making the 
parameter renormalization) and the lower end of the physical flow time t is limited by the 
lattice spacing as a <C \/^. Thus the asymptotic t —>• 0 behavior in Eq. (ini) will be useful 
to find the extrapolation for t —>• 0. 


^ For completeness, we quote the related formulas: The running gauge coupling is defined by 


^ dg(g) 
^ dg 


/3(5(9)), 


ff(9 = g)=9, 


( 1 . 8 ) 


where the beta function is given from the renormalization constant Z in the MS scheme in 

€6^ 
by 


P 2 ^ ( 1 ) 

-=-^-9 ■ 

9 d,9 


(1.9) 

( 1 . 10 ) 
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The above universal formula for the axial-vector current is quite analogous to universal 
formulas for the energy-momentum tensor—the Noether current associated with the trans¬ 
lational invariance—on the basis of the gradient/Wilson flow 16Hl9l|: see also Ref. 20(]. 
The motivation is also similar: Since lattice regularization breaks the translational invari¬ 
ance, the construction of the energy-momentum tensor with lattice regularization is not 
straightforward 2l|, [2^. The universal formulas in Refs. lfiHl9l|. because they are thought 


to be regularization independent, may be used with lattice regularization. The validity of 
the formulas has been tested W employing Monte Carlo simulations 23|, [2^ and the 1/N 


expansion 


25( (see also. Ref. 2^ for the 1/N expansion of the gradient flow). 

A main technical difference in the derivation of Eg. (11.71) f rom the derivation of the universal 
formulas for the energy-momentum tensor in Refs. |l6l4l9l] is that only the knowledge of the 
one-loop expression in dimensional regularization (jl.5|) will be used in what follows. On the 
other hand, in Refs. lfiHl9l|. full-order perturbative expressions of the correctly normalized 
energy-momentum tensor (which is readily obtained by dimensional regularization) were 
used. Although it must be possible to arrive at Eq. dm) starting from an “ideal” axial- 
vector current obtained from the chiral-symmetry-preserving lattice fermions by the Noether 
method, a much simpler one-loop egression such as Eq. (jl.5p is sufficient. This point, we 
think is technically very interesting^ 

2. Gradient/Wilson flow 

The gradient/Wilson flow [3-0] is an evolution of field configurations according to flow 
equations with a fictitious time t. Eor the gauge field A^(x), the flow is defined by [3, 0] 

dtBfj,{t,x) = D„Gufj.{t,x) + aoD^d^B^{t,x), Bf,{t = 0,x) = A^{x), (2.1) 

where Gfj,u{t,x) is the field strength of the flowed gauge field R^(t,x), 

G^y{t,x) = dfj_B„{t,x) - duBf,{t,x) + [Bf,{t,x), B^{t,x)], (2.2) 

and the covariant derivative on the gauge field is 

D/j, = d/j, + [B/j,, •]. (2-3) 

The second term in the right-hand side of Eq. (j2.1|) is a “gauge-fixing term” that is introduced 
to provide the Gaussian damping factor (see below) to the gauge mode. Although it breaks 
the gauge covariance, it can be shown that any gauge-invariant quantity is independent of 
the “gauge parameter” ao |3|- 
Eor the quark fields, '<p{x) and 'ip{x), the flow is defined by 

dtx{t, x) = [A- aodf,Bf,{t, x)] x{t, x), x(t = 0,x) = ^(x), 


dtx(t,x) = x{t,x) "a +aodi,Bi,{t,x) , x{i = 0,x) =-i/ix), 
where covariant derivatives for the flowed quark field, xi^^x) and x{t,x), are 

A = D^D^, = d^ + Bfj^, 


(2.4) 

(2.5) 

( 2 . 6 ) 
(2.7) 


^ The idea that the knowledge of a one-loop expression such as Eq. (11.51) would be sufficient to find 
the universal formula has emerged through discussion with David B. Kaplan. We would like to thank 
him. 
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For the implementation of these flow equations in lattice gauge theory, see Refs. 


and are quantum fields being 
perturbation theory for quantum 


The initial values in above flow equations, A^{x), il:[x 
subject to the functional integral. One can develop 
correlation functions of the flowed fields, x), x(t, x), and x(t, x). For example, the tree- 
level propagator of the flowed gauge field (in the “Feynman gauge” in which Aq = ao = 1, 
where Aq is the conventional gauge-fixing parameter) i.'^ 


J2 srab s 


J V 


Similarly, the tree-level quark propagator is 

{X{t,x)x{s,y)) = [ 

Jv 


ip{x-y) 


=-{t+s)P^ 

p2 


-{t+s)p‘^ 

ij) -I- Mo 


(2.9) 


( 2 . 10 ) 


The details of the perturbation theory (the flow Feynman rule) are summarized in Ref. 171. 

A remarkable feature of the gradient/Wilson flow is its ultraviolet (UV) finiteness [3|, l4|. 
Correlation functions of the flowed gauge field become UV finite without the wave function 
renormalization (if one makes the parameter renormalization). This finiteness persists even 
for local products, i.e., composite operators. The flowed quark field, on the other hand 
requires the wave function renormalization. However, once the elementary flowed quark field 
is multiplicatively renormalized, any composite operator of it becomes UV hnite. Thus, if 


one introduces the combinations 


Szl 


X{t,x) 


X{t,x) 


—2 dim{R)Nf 


{^Xit,x) px{t,x] 


—2 dim(i?)Vj 




^ {47r)H‘^ (^X{t,x) px{t,x) 


X(t,x), 


Xit,x), 


( 2 . 11 ) 

( 2 . 12 ) 


where 

"0^ = 0^- (2.13) 

then wave function renormalization factors are canceled out in x(i) x) and x(ij x) and any 
local product of them becomes a UV-finite renormalized operator. Thus, if one can express 
a composite operator in terms of local products of flowed fields, B^{t, x), x(t, x), and x(i, x), 
then the expression provides a universal formula for the composite operator. What we shall 
do is to apply this idea to the composite operator in Eq. m- An explicit method to rewrite 
a composite operator in terms of local products of the flowed fields is given by the short 
flow-time expansion that is the subject of the next section. 


’ Throughout this paper, we use the abbreviation 




(27r) 


D ■ 


( 2 . 8 ) 


^ Here dim(i?) is the dimension of the representation R of the gauge group to which the fermion 
belongs (dim(i?) = 3 for quarks). The coefficients in these expressions are chosen so that x{t,x) —>■ 
X(t,x) and x{t,x) x(bx) for < 0 [nl- 
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3. Small flow-time expansion 

We take a would-be axial-vector current composed of the flowed quark field 

(3.1) 

and consider its t —)■ 0 limit. As discussed in Ref. j^, a local product of flowed fields in the 
t —)• 0 limit can be expressed by an asymptotic series of local composite operators of fields at 
zero flow time with increasing mass dimensions. In the present case, because of symmetry, 
we have 

X{t, x)jf,j5t^x{t, x) = c{t)il;{x)jf,X5t^ip{x) + 0{t) (3.2) 

for t —)■ 0. 

The expansion can be worked out in perturbation theory; the coefficient c{t) can be found 
by computing the correlation function 

(x(i, x)'ilj{y)'4>{z)) (3.3) 

and comparing its t —0 behavior with 

{'4>{x)xi,X5t'^^p{x)^p{y)^p{z)) . (3.4) 

To one-loop order, there are seven IPI “flow Feynman diagrams” which contribute 
to Eq. (|3.3p : they are depicted in Figs [THZ] according to the convention in Ref. [I^. Among 


Fig. 1 Diagram a 





these diagrams, diagrams b, c, and e in Figs. [5HZ1 are irrelevant for the coefficient c(t) 
in Eq. (13.21) because they are proportional to the momentum of the external quark line. 
Then, writing Eq. (13.31) as 



d^q 

{2-k)D 


^ip{y-x)^iq(x-z) 




(3.5) 


the contribution of diagram a is A4^^(0,0; t) = For one-loop diagrams, we work with 

dimensional regularization [with the prescription (11.21) ]. The contribution of each diagram 
to A4^^(0,0;t) is tabulated in Tabled) Combining all the contributions, we have 
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Fig. 5 Diagram b 


Fig. 6 Diagram c 


Fig. 7 Diagram e 


Table 1 Contribution of each diagram to 0;t) in Eq. (13.511 in units 


d 

f 

g 


(- 1 ) 


1 7 

- + ln(87rt) + - 
€. A 


1 


+ ln(87rt) + 1 


(-4) 


1 , . N 1 

- + ln(87rt) + - 
6 ^ 


.,2 

(dvr)^ 


A4^„(0,0; t) = <1 1 + -^C2(fi)(-3) - + ln(87rt) + \ 


6 


7/^75^^ 


(3.6) 


Invoking a facile method explaining in Ref. 271], this implies the small flow-time expansion 
is given by 


-.2 

(47r)2 


= <1 1 + 7^C'2(R)(-3) i+ln(87rt) + J 


+ o{t). 


(3.7) 


Next we express this in terms of the “ringed” fields in Eqs. (12.1111 and (12.1211 . Using the 
result of Ref. [l3| that 


—2dim(R)A^f 

{ATr)H‘^ (^X{t,x)'^x{t,x] 

where 


= Zie) 1 + 


9 


(dvr)' 


rC'2(R) 


3- - ^{t) 
e 


+ 0{g^) -I- 0{t) > , (3.8) 


Z{e)^ 


{SirtY ’ 


^{t) = —3 ln(87r//^t) -|- ln(432), 


(3.9) 


we have for e —)• 0, 


x{t,xhi,X5t^x{t,x) 


= u + 


9 


tC2{R) 


(47r)2 

where we have set gQ = 


---ln(432) 


+ 0(5 ) } Y{x)x^ilbt Y{x) + 0(t), 


(3.10) 
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Inverting this relation for t —)• 0 and plugging it into Eq. (HH), we have 




--+ln(432) 


+ 0(5 ) } x{t, a:)7M75i x(i, x) + 0{t). (3.11) 


As the Ward-Takahashi relation (jl.ll) shows, j^^{x) in the left-hand side of this relation 
must be UV finite. The right-hand side of this relation is certainly UV finite, as all 1/e 
singularities are cancelled out. 

Finally, we invoke a renormalization group argument. By applying the operation 



(3.12) 


to the both sides of Eq. dni), where the subscript 0 implies the bare quantities are kept 
fixed, we have (/x9/0/r)oj^(x) = o|§We also have {iJ,d/dfj,)ox{t,x)jf^'y 5 t^x{t^x) = 0, because 
the flowed quark fields are certain (although very complicated) combinations of bare fields. 
Therefore, the quantity in the curly brackets in Eq. (13.1011 is independent of the renormal¬ 
ization scale /r. This implies that if one uses the running gauge coupling g{q) in place of g, 
the expression is independent of the scale q. Thus, we set q = 1/V^ as a particular choice. 
Since g{l/V^) —^ 0 for f —>■ 0, the perturbative computation is justified in the f —)■ 0 limit 
which also eliminates the 0{t) term in Eq. (|3.1ip . In this way, we arrive at the universal 
formula dlZD. 

An argument similar to the above may be repeated for the flavor non-singlet pseudo-scalar 
density that fulfills the PCAC relation. 


= -<5^(y - x)'y 5 t^ {ip{y)ip{z)) - 5‘^{z - x) {ip{y)'4}{z)) 75 ^^. (3.13) 


Omitting all details of the calculation, the final result is given by 

^{x) 

^ + ^^^(4^^ C* 2 (-R) [31n7r-b6-Fln(432)]| x(i,x) 75 {t^,M(l/v^)}x(t,a:). 

(3.14) 


In these expressions, the MS or MS scheme is assumed both for the renormalized mass 
matrix M and the renormalized pseudo-scalar density; M(q) is the running mass matrix at 
the renormalization scale q. 


4. Triangle anomaly 

The matrix element of the above flavor non-singlet axial-vector current is relevant for the 
pseudo-scalar meson decay process by weak interaction. For the process vr^ —)■ 27 , the triangle 
diagram containing one axial-vector current and two vector currents is important 2^, |29(] . 


® Since j^^{x) is the unique gauge invariant flavor non-singlet dimension 3 axial-vector operator, 
{gd/dg)od^j^j^{x) = 0 implies this. 
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It is thus of great interest to examine the three-point functioi|§ 


■ABC 


T 

'fiup 


{p,q) = I d^x I d^ye' 


^D..Jpx iqy 


7 ix)jf {y)) + iB^ C)] 


(4.1) 


where i^(0) is defined by our universal formula (11.71) and and j^{y) are flavor non¬ 

singlet vector currents 


jv (y) = i’{y), jp {z) = 'tPizhpt ipiz). 


(4.2) 


Usually, conventional regularization preserves the flavor vector symmetry and a naive 
expression of the vector current is correctly normalized. 

In the present paper, we investigate whether (|4.ip reproduces the correct triangle (or axial) 
anomaly in the lowest non-trivial order of perturbation theory. We work with massless theory 
for simplicity. In the lowest order of perturbation theory, using Eqs. (11.711 . (|2.1ip . (I2.12p . (|3.8p . 
and hnally (j2.inp . Eq. (|4.ip before taking the t —)• 0 limit is given bjo 


7;^p^(p,g)|o(g0) = -7trij(l)tr 


X / tr 'Y o-t; ^ ^ o-t{(-+pYp-d^-qY 

X / tr^7^75^^^7.^7p^_^e e 


1 1 1 






(4.3) 


The total divergence of the axial-vector current, i.e., the triangle anomaly, is thus given by 


i{P+(l)p 'Vp^(P,9)|o(gO 


= -^tr/j(l)tr 


X / trys 

Je 


1 


1 




p). 


(4.4) 


It is interesting to note that if there were no Gaussian damping factors which result from 
the flow of the quark fields in this expressioii, a naive shift of the loop momentum makes 
this expression vanish as the well-known case 28|, l29( . The reality is that there are Gaussian 
damping factors and, in the f —)• 0 limit, we have the following non-zero result: 


p —2tP 

*(P + 9)M'7)i^P^(p,9)|o(gO) = -2t y^^-^trR(l)tr ,t^}) e^uppPaq^ 

= ^auppPaqp- (4.5) 

However, this is not quite identical to the conventional triangle anomaly; the coefficient is 
half the conventional one. We recall that the coefficient of the conventional triangle anomaly 


®Here, we supposed that = 0. 

^°tri{(l) denotes the trace of the unit matrix over the gauge representation index; tr/{(l) = 3 for 
quarks. 
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is fixed by imposing the conservation law of the vector currents. Being consistent with this 
fact, we observe that the vector current is not conserved with Eq. (14.,'ll) : 

- ipy (^afi/3pPaqi3- (4.6) 

Do these observations imply that our universal formula (1121) is incompatible with the 
physical requirement that the vector currents are conserved (i.e., vector gauge invariance)? 

We should note, however, that it is not a priori clear whether the small flow-time expan¬ 
sion (13.2p for t —>• 0 holds even if the point x collides with other composite operators in 
position space (see the discussion in Sect. 4.1 of Ref. [l^, for example). The integration (14.11) 
in fact contains the correlation function at equal points, a; = 0 or y = 0. Therefore, there 
exists freedom to modify Eq. (BSD by adding a term that contributes only when x = 0 
or y = 0 in position space. Using this freedom, we can redehne the correlation function as 

{j 5 pi^)ju ix)jp (y)) - Y^trij(l)tr [dj'^{x)5'^{y) - 6^{x)da6{y)] . 

(4.7) 

This redefinition preserves the Bose symmetry among the vector currents and, of course, 
does not affect the correlation function when x 7 ^ 0 or y 7 ^ 0 . 

After the redefinition, we have 

iiP + q)p '7p^p'^(P,^)|o(gO) = -^trij(l)tr eauPpP^qp, (4.8) 

- iPu 7))tp^(P,g)|o(gO) = -iqp 7p^p^(P,g)|o(gO) = 0. (4.9) 

These expressions coincide with the conventional form of the triangle anomaly. Since what 
we added in Eq. ( 021 ) is simply a term that vanishes for x 7 ^ 0 or y 7 ^ 0 , our computation 
above shows that the universal formula (II.7p produces non-local structure that is consistent 
with the triangle anomaly, at least in the lowest non-trivial order of perturbation theory. 


5. Conclusion 

In this paper, we presented a universal formula that expresses a correctly normalized flavor 
non-singlet axial-vector current through the gradient/Wilson flow. The formula is universal 
in the sense that it holds irrespective of the chosen regularization, and especially holds with 
lattice regularization. Whether our formula possesses possible advantages over past methods 
in actual lattice Monte Carlo simulations must still be carefully examined. In particular, 
the small t limit in Eq. ( 021 ) is limited by the lattice spacing as a <C and there exists a 


systematic error associated with the extrapolation for t —)• 0. See Refs. [23|] and 2^ for the 
t —>• 0 extrapolation in Monte Carlo simulations with the universal formula for the energy- 
momentum tensor. 

As a purely theoretical aspect of our formula, it is interesting to note that if one can show 
that the triangle anomaly obtained by i{p + q)p.Ti^p^{p,q) from Eq. (14.ip is local in the 
sense that it is a. polynomial of p and q (that we believe is quite possible), then we may 
repeat the proof 30[ of the Adler-Bardeen theorem i.e., the triangle anomaly does not 
receive any correction by strong interaction. 
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It must also be interesting to generalize the construction in the present paper to the flavor 
singlet axial-vector current. Here, one has to incorporate the mixing with the topological 


charge density. Then it is quite conceivable that this construction give a further insight on 
the nature of the topological susceptibility dehned through the gradient/Wilson flow [21. l32|. 


Acknowledgments 

Participation at the KEK lattice gauge theory school organized by JICFuS was quite helpful 
to us and we would like to thank the organizers for their hospitality. The work of H. S. is 
supported in part by Grant-in-Aid for Scientific Research 23540330. 

References 

[1] S. Aoki, Y. Aoki, C. Bernard, T. Blum, G. Colangelo, M. Della Morte, S. Diirr and A. X. El Khadra 
et ai, Eur. Phys. J. C 74, no. 9, 2890 (2014) arXiv: 1310.8555 [liep-lat]]. 

[2] M. Liischer, JHEP 1008, 071 (2010) [Erratum-ibid. 1403, 092 (2014)] arXiv:1006.4518 [hep-lat]]. 

[3] M. Liischer and P. Weisz, JHEP 1102, 051 (2011) arXiv:1101.0963 [hep-th]]. 

[4] M. Liischer, JHEP 1304, 123 (2013) arXiv: 1302.5246 [hep-lat]]. 

]5] A. Shindler, Nncl. Phys. B 881, 71 (2014) arXiv: 1312.4908 [hep-lat]]. 

]6] D. B. Kaplan, Phys. Lett. B 288, 342 (1992) hep-lat/9206013 . 

]7] Y. Shamir, Nucl. Phys. B 406, 90 (1993) hep-lat/9303005 . 

]8] V. Furman and Y. Shamir, Nncl. Phys. B 439, 54 (1995) hep-lat/9405004 . 

[9] H. Neuberger, Phys. Lett. B 417, 141 (1998) hep-lat/9707022 . 

[10] H. Neuberger, Phys. Lett. B 427, 353 (1998) hep-lat/9801031 . 

[11] P. H. Ginsparg and K. G. Wilson, Phys. Rev. D 25, 2649 (1982). 

[12] P. Hasenfratz, Nucl. Phys. B 525, 401 (1998) hep-lat/9802007 . 

[13] P. Hasenfratz, V. Laliena and F. Niedermayer, Phys. Lett. B 427, 125 (1998) hep-lat/9801021 . 

[14] M. Liischer, Phys. Lett. B 428, 342 (1998) hep-lat/9802011 . 

[15] J. C. Collins, Renormalization. An Introduction to Renormalization, the Renormalization Group, and 
the Operator Product Expansion (Cambridge University Press, Cambridge, UK, 1984). 

[16] H. Suzuki, PTEP 2013, no. 8, 083B03 (2013) arXiv: 1304.0533 [hep-lat]]. 

[17] H. Makino and H. Suzuki, PTEP 2014, no. 6, 063B02 (2014) arXiv:1403.4772 [hep-lat]]. 

[18] H. Makino and H. Suzuki, arXiv: 1410.7538 [hep-lat]. 

[19] H. Suzuki, arXiv:1501.04371 [hep-lat]. 

[20] L. Del Debbio, A. Patella and A. Rago, JHEP 1311, 212 (2013) arXiv:1306.1173 ]hep-th]]. 

[21] S. Caracciolo, G. Curci, P. Menotti and A. Pelissetto, Nucl. Phys. B 309, 612 (1988). 

[22] S. Caracciolo, G. Curci, P. Menotti and A. Pelissetto, Annals Phys. 197, 119 (1990). 

[23] M. Asakawa et al. [FlowQCD Collaboration], Phys. Rev. D 90, no. 1, 011501 (2014) arXiv:1312.7492 
[hep-lat]]. 

[24] M. Kitazawa, M. Asakawa, T. Hatsuda, T. Iritani, E. Itou and H. Suzuki, arXiv: 1412.4508 [hep-lat]. 

[25] H. Makino, F. Sugino and H. Suzuki, arXiv: 1412.8218 [hep-lat]. 

[26] S. Aoki, K. Kikuchi and T. Onogi, arXiv: 1412.8249 ]hep-th]. 

[27] H. Makino and H. Suzuki, Appendix D of arXiv: 1403.47721^4 [hep-lat]. 

[28] S. L. Adler, Phys. Rev. 177, 2426 (1969). 

[29] J. S. Bell and R. Jackiw, Nuovo Cim. A 60, 47 (1969). 

[30] A. Zee, Phys. Rev. Lett. 29, 1198 (1972). 

[31] S. L. Adler and W. A. Bardeen, Phys. Rev. 182, 1517 (1969). 

[32] J. Noaki et al. [JLQCD Collaboration], PoS LATTICE 2013, 263 (2014). 


11 





